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Abstract 

We shall study a turbulence model arising in compressible fluid mechanics. 
The model called 9 — cp we study is closely related to the k — e model. We 
shall establish existence, positivity and regularity results in a very general 
framework. 
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1 Introduction 

We shall first recall some basic ideas concerning the modelisation of the tur- 
bulent fluids, the reader can consult [23l|27j for a more detained introduction. 



Let u,p,p,T be the velocity, pressure, density and temperature of a new- 
tonien compressible fluid. Let also !1cl 3 a domain which is asummed to be 
bounded. Then the motion of the flow in O at a time t 6 M + can be described 
by the compressible Navier Stokes equations (see system (C) page 8 in [18]). 
It is well known that direct simulation based on such a model is harder or 
even impossible at high reynolds numbers. The reason is that too many points 
of discretization are necessary and so only very simple configurations can be 
handled. 



Thus engineers and physicists have proposed new sets of equations to de- 
scribe the average of a turbulent flow. The most famous one is the k—e model, 
introduced by Kolmogorov [13] . We shall briefly present its basic principles in 
the following. Let v denote a generic physical quantity subject to turbulent 
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(i.e. unpredictable at the macroscopic scale), we introduce its mean part (or 
its esperance) (v) by setting: 

(v(x,t)) = / v(oj,x,t)dV(ui), 
Jv 

where the integral is taken in a probalistic context which we shall not detail 
any more here. Note however that the operation (.) is more generally called a 
filter. The probalistic meaning is one but not the only possible filter (see for 
instance |23| chap. 3). We shall then consider the decomposition: v = (v) +v' , 
where v' is refered to the non computable or the non relevant part and {v } is 
called the mean part (i.e. the macroscopic part). 

The principle of the k — e model is to describe the mean flow in terms 
of the mean quantities (u) , (p) , (p) , (T) together with two scalar functions k 
and e, which contains relevant informations about the small scales processes 
(or the turbulent processes). The variable k (SI: [t-]) is called the turbulent 
kinetic energy and e px] is the rate of dissipation of the kinetic energy. They 
are defined by: 

k = l^\ 2 £ = |<|Vu' + (Vu') T | 2 >, (1) 

where v is the molecular viscosity of the fluid. The model is then constructed 
by averaging (i.e. by appling the operator (.) on) the Navier-Stokes equations. 
Under appriopriate assumptions (i.e. the Reynolds hypothesis in the incom- 
pressible case, and the Favre average in the compressible case) we obtain a 
closed system of equations for the variables (u), (p), (p), (T), k and e (see [23] 
pages 61-62 for the incompressible case, and pages 116-117 in the compressible 
situation). 

Here we shall focus on the equations for k and e and we consider that the 
others quantities (u), (p), (p), (T) are known. Moreover, in order to simplify 
the readability we do not use the notation (.), i.e. in the sequel we will write 
u instead of (u) and p instead (p) to represent the mean velocity and density 
of the fluid. The equations for k and e are of convection-diffusion-reaction 
type: 

dtk + u-Vk- — div {(u + p—)Vk) = c u —F - -kD - e, (2) 

p e e 3 

c A: 2 2ci e 2 

d t e + u • Ve - —div (h> + p— )Ve) = c x kF -eD - c 2 — , (3) 

p e 6c v k 

where D(x,t) := divu(x,t), F(x,t) := ^|Vu+ (Vu) T | 2 - §L>(:r,t) 2 > (see 
subsection l6.2l in the Appendix) and c u ,c £ , c\,C2 are generally taken as positive 
constants (see (|70]) in the Appendix). 

Note that equations d2])- (J3j) are only valid sufficiently far from the walls. 
In fact, in the vicinity of the walls of the domain 0, there is a thin domain E, 
called logarithmic layer in which the modulus of the velocity goes from to 
0(1). In this layer we can use some wall law or a one equation model (see [23J 
chap.l and [24]) instead of (ED- (J3j) . Note however that the equations (EJ) — (J3J) 
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can be considered even in the logarithmic layer if we allow the coefficients 
c u ,c £ ,ci and C2 to depend appriopriately on some local Reynolds numbers 
(see [23J pages 59-60 and page 115). In this last situation the system is called 
Low-Reynolds number k-e model. 

In the following we focus on the study in the domain f2 := f2 \ S and we 
assume that its boundary dQ is Lipschitqll. We denote by n(x) the outward 
normal defined for all points x G d£l. The boundary conditions for k and e 
on d£l are well understood. We have: 

k = ko and £ = £o on dQ, (4) 

where ko and £o are strictly positive functions which can be calculated by 
using a wall law (see [23] p. 59) or a one equation model (see [21] )• In the 
following we assume that ko and £o are given. Moreover we can assume (see 
again [23] p. 59) that: 

u • n = on d£l. (5) 

We shall concentrate in this paper on a reduced system called 9 — (p model. 
The new variables 9 [s] and if [m~ 2 ] are obtained from k and e by the formulas: 

k F 2 

These variables have a physical meaning (see [27J): 9 represents a characteris- 
tic time of turbulence and L = y? -1 / 2 is a characteristic turbulent length scale. 
By using this change of variable in the equations d2J)- ([3]) and after considering 
some modelisation arguments for the diffusion processes (see the Annexe) we 
obtain: 



(P) 



where the coefficients cq,c v and Cj are all positive. 

Problem (P) is known as the 9 — ip model. It differs from the k — e one 
only by the diffusive parts and it is attractive by some stronger mathematical 
properties. Another model closely related to these systems, and having some 
popularity, is the k — ui one (with uj = 9~ l , see for instance [27]). 
In the papers [IH [20] the authors have established the existence of a weak 
solution for (P) and a property of positiveness. This last feature takes the 
model useful in practice: it can be used directly or also as an intermediary 
stabilization procedure to the k — e one (see [22]). Another important prop- 
erty attempted for a turbulence model is its capability to predict the possible 
steady states. In the previous works only the evolutive version of (P) was 
studied (except in [21] where however only the incompressible situation, with 
pertubated viscosities was considered), and the results obtained cannot pre- 
dict the existence or non-existence of steady states. 

Hence in this paper we shall study the stationary version of (P) on a 
bounded domain C Mr, N=2 or 3, on which we impose the boundary 

x see [9] p. 127 



d t 6 + u • V0 - i div ({v + if)Wj = -c 3 6 2 F + Ci 9D + c 5 

dt(p + u • Vip - \ div Uy + ^)Vv?) = - V (c 6 9F - c 7 D + cg^ 1 ) 
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conditions 9 = a, p = b on d£l. Remark that by using @ together with © 
we obtain 

a = T b = ¥ (7) 

Hence we can assume that o and b are strictly positive given functions. 

We shall establish existence, positivity and regularity results in a very 
general framework. In [8] we established existence and regularity results for a 
turbulent circulation model involving u and k as unknowns. The reader inter- 
ested for recent references concerning the numerical simulation of turbulent 
fluid can consult 25 1. 



2 Main results 

2.1 Assumptions and notations 

Let (Q) denote the stationary system associated to (P). For simplicity we in- 
troduce the new parameters C\ n & := pc{ n d where the subscript 'ind' takes the 
integer values 3,4,5,6,7,8 or the letters 9 and (p. Then our main model (Q) 
has the following form: 



(Q) 



pu-VO- div ({y + )V0) = -C 3 F9 2 + C A 9D + C 5 in Q 



6tp< 

pu-Vip- div ( (y + ^)Vv? ) = -(p{Ce6F - C 7 D + C^ -1 ) in n 



a, if = b on 8ft 



For physical reasons we are only interested in positive solutions (9, ip) for (Q). 
Note however that even with this rectriction, the problem (Q) may be singu- 
lar (i.e. the viscosities v + and v + may be unbounded). Moreover, 
because we allow v = the equations may degenerate (i.e. the viscosities 
may vanish). Hence without additional restriction there may be various non 
equivalent notions of weak solution (see for instance [8]). 

In fact a good compromise between respect of the physics, simplification 
of the mathematical study and obtention of significative results, is to restrict 
9 and ip to be within the classe S defined by: 

S = {/ : Q R+ such that / G H 1 ^) n L°°(fi), / _1 E L°°(ft)}. 

In particular, if the parameters appearing in (Q) are sufficiently regular and 
if we restrict 9 and p to be within the classe S, then the notion of a weak 
solution for (Q) is univocally defined: it is a distributional solution (9, ip) that 
satisfies the boundary conditions in the sense of the trace. 
In this last situation we will tell that (9, ip) is a weak solution of (Q) in 
the class S. 

In order to can consider such a weak solution for (Q) we shall precise in 
the following some sufficient conditions of regularity for the data. 
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Let N=2 or 3 denote the dimension of the domain Q, and r be a fixed number 
such that: 

N 

r > y (8) 
We then have the following continuous injection (see lemma E]): 

L r (Q) C W~ l ' p (p), where (3 = r* > N. (9) 

Recall that D = div(u). We will consider the following assumptions: 

• Assumptions on S7: 

0, C 1^ is open, bounded and it has a Lipschitz boundary dQ. (10) 

• Assumptions on the flow dataElu, F,D,p and u: 

i/>0, (11) 

F,p: n^R + , p, p^ 1 G L°°(p,), F G L r (Q), (12) 

u G (L 3 ^))^, £> G L r (0), div(pu) = 0, (13) 

u • n = on dn. (14) 

• Assumptions on the turbulent quantities on the boundary: 

a,be H 1/2 {dn)C\L°°{dti), 

a(x), b{x) > 5 > a.e. x G d$l, (15) 

where S > is a fixed number. 

• Assumptions on the model coefficients: 

d : O R+ C 5 , C 8 G L r (^), for t / 5, 8 : d G L°°(ft) (16) 

C ind :Ox(I + f^K + are Caratheodory (17) 

Cind(x,v,w) > ai n d > \/v,w G M + and for a.a. (18) 

Vu, w G M + , x -> C ind (x, u, w) G L°°(n), (19) 

where in (ri~7]) - (fT9j) Ci n d means or C v . The assumption (JTTj) signifies 
that \/v,w G M + , x — > C- m ^{x,u,v) is measurable and for a.a. x G : 
(v,?if) — )• Ci n d(x, f , to) is continuous. This ensures that Ci n d{x , 9 , ip) is 
measurable when 9 and </? are measurable. The condition ()18f) means that 
Ci n d is uniformly positive, whereas (fl~9j) tells that Ci n d(x, 6*, </?) remains 
bounded if 9 and 99 are bounded. 

We will study problem (Q) under the main assumption: 

(Aq) : © - (USD are satisfied. 

Note that in the main situation (Aq) the assumption (jlip made for v al- 
lows the possibilty v = 0. In other words the molecular viscosity v can be 
neglected in the model. This is often chosen in practice because the eddy 
viscosities and are dominant in the physical situations (see 



2 when N = 2 one assumption in (fl3|) can be relaxed: u e (L p (il)) 2 with p > 2 (instead of 
p = 3) is sufficient, but this would not improve any result significantly. 



5 



Remark also that the coefficients C, are allowed to depend on x, and the vis- 
cosity parameters Cg, C v may depend on x,0,ip. 

For a given function / : 0, — > R, we shall use the notations / + and /~ to 
represent the positive and negative parts of /, that is: 

f = f + + r, f + (x)=m a x(f(x),0)>0, f-(x) =min(/(x),0) < 0. 

We will also consider some assumption of low compressibility of the form: 

\\D + \\ L r(n) < t, (20) 

for some r > that will be precised. 

This last kind of condition seems to be necessary (see the Appendix) in order 
to obtain a weak solution for (Q) in the three dimensional case, whereas when 
N = 2 we shall use some particularities of the situation to obtain a weak 
solution without any assumption of low compressibility. Nevertheless in this 
case we will assume that in addition to (Aq) the following condition is satisfied: 

C 2 

(Ax) : -f G L r , c 6 = 0.014, c 7 = 0.104, c 8 = 0.84. 

^3 

In this last condition the values for cq , cj and cq are in fact their classical 
constant values (see ([75]) in the Appendix) 

In the sequel we denote by DATA some quantity depending only on the 
data under the assumption (Aq), i.e.: 

DATA = Constat, a, b, \\u\\^N,ae,a v , (\\Cj\\ L r)j =5>8 , \\F\\ L r, 

(||Ci||L 00 )i=3,4,6,7)- (21) 

Note that DATA does not depends on D and v. 

The exact form of the dependency (i.e. the function Const) is allowed to 
change from one part of the text to another. 

2.2 Main results 

We shall establish two theorems of existence. The first one applies if N=2 or 
3 and the second one is limited to the case N=2. 

Theorem 1 Assume that (Aq) holds. Then there exists r > such that if 
||D + ||^r(Q) < t then problem (Q) admits at least one weak solution (8,ip) in 
the class S. 

Theorem 2 Assume that N = 2 and (Aq),(A\) hold. Then problem (Q) 
admits at least one weak solution (9,ip) in the class S. 

In all the situations we have the following regularity result: 

Theorem 3 Let (9, ip) be a weak solution of (Q) in the class S and assume 
that (Aq) is satisfied. We have: 

i) Ifu G (L°°(Q)) N anda,b are Holder continuous, then(9,cp) G (C°' a (H)) 2 , 
for some a > 
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ii) Assume that in addition, the following conditions are satisfied: 

dQ, a and b are of class C 2 ' a , pu G (C 1 '"^))^ 
F G C°> a {n), \/i = 3, .., 8 : Q G C > a (Ti) 
C md (x,v,w)eC 1 > a (nx(R + f) 

Then (0,ip) G (C 2,a (f2)) 2 and it is a classical solution of (Q). 

2.3 Discussion on the results 

The major feature of our work is to treat the compressible situation in a gen- 
eral framework. The compressible case is interesting for several applications 
(see for instance [22L [5]). However the additional terms of the model arising 
from the compressibility effects induce important complications for its analy- 
sis. Roughly speaking: these additional terms are responsible of an increase 
or decrease (depending on the signe of D) of the turbulence. Then the bal- 
ance between the increase/decrease of the source terms of turbulence and the 
possible explosion/ vanishing of the turbulent viscosities is difficult to control. 
Compared to previous works (see [HI [151 1231 EE]) our basic assumption (Aq) 
made in theorem [T] is very general, and we remark that: 

i) We do not impose a free divergence condition on u and the mean density 
p is not supposed constant. Hence our analysis can handle compressible 
turbulent flows. The condition of the form (j20f) assumed for divu is 
much weaker than the corresponding assumption made in [15]. In di- 
mension two, under the additional condition (^4i) we obtain in theorem 
E] an existence result without any additional assumption of low com- 
pressibility. The necessity of a condition of the form (I20jl when N=3 is 
discuted in the Appendix. 

ii) We allow the viscosity parameters Cq^C^ to depend on x,8,ip whereas 
in the previous works these parameters were taken constant (except in 
|20[ [2~Tj where however an artificial regularization was made, and only 
the incompressible situation was studied). 

iii) We only assume weak regularity for the mean flow, i.e: u G (L^(tl)) N 
and divu G L r (£l), with some r > 3/2, whereas in the previous works 
it was assumed u G (L°°(Q)) N and divu G L°°(f2). Our assumption is 
more interesting from a mathematical point of view because it is satisfied 
when u is a weak solution of the Navier Stokes equations. Hence our 
work could be used for a future analysis of the full coupled system Navier- 
Stokes plus (Q). 

iv) All the coefficients Ci are allowed to depend on x G £1. Hence our study 
is also a step for the analysis of the Low- Reynolds k — e model (in [28] 
the Low-Reynolds k — e model is studied only in the incompressible and 
unsteady situation). Finally the boundary values for 9 and (p are not 
taken constants unlike the previous works. 

Note also that under the additional assumption u G (L°°(Q)) N we give 
a Holder continuity result for k and 8. Moreover, we establish an existence 
result for a classical solution under some smoothness assumptions on the data. 
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These regularity results seem to be completly new. 



Another feature of our work is to point out that the choice (p 




indicated even in the compressible situation: this makes the dynamic of the 
(p equation stable when N = 2 and stable under additional conditions when 
N = 3 (see the Appendix). These results improve the study of the model 
presented in [23J, chap. 9. 

In order to establish our theorems we establish intermediate results (see 
Proposition [9]), concerning the existence, positivity and regularity properties 
for a weak solution of an elliptic scalar problem (possibly degenerate and 
singular) of the form: 



These results have also an independent mathematical interest. 

2.4 Organization of the paper 

• In section 3 we shall recall some results concerning: the truncature at a fixed 
level and the Stampacchia's estimates. This last technique takes an important 
role in our analysis, moreover we shall need a precise control of the estimates. 
Hence we shall present it with some details and developments. 

• In section 4 we introduce a sequence (Q n ) of problems which approximate 
(Q). For n fixed (Q n ) is a PDE system of two scalar equations of the form (S): 
one equation for the unknown 9 n +i and one for (p n +i- The point is that the 
unknowns 9 n+ \ and fn+i are only weakly coupled. The coupling of the two 
equations is essentially realized through the quantities 9 n and ip n calculated at 
the previous step. Hence we shall firstly study carefully the problem (S). The 
major tool used here are the Stampacchia's estimates. We next prove that 
(Q n ) is well posed. Hence we obtain an approximate sequence of solutions 
(8 n ,ip n ) for problem (Q). Moreover, we prove that 9 n and ip n are uniformly 
bounded from above and below, which are the key estimates. 

• In section 5 we use the uniform bounds established in section 4, in order to 
extract a converging subsequence from (9 n , <p n ). We then prove that the limit 
(9, ip) is a weak solution of (Q) in the class S. 

Under the additional assumption u 6 (L°°(f2)) 3 we are able to use the De 
Giorgi-Nash Theorem and we obtain an Holder continuity result for 9, (p. By 
assuming in addition some smothness properties for the data we can iterate 
the Schauder estimates and prove Theorem O 

• In section 6 (Appendix) we present the derivation of the 9 — ip model from 

2 

the k — e one. Moreover we justify that the choice ip = is valid even in 
the compressible situation. The justification uses in particular a property 
of positivity of the function F. We also dicuss briefly the necessity of the 
low compressibility assumption when N = 3. Finally we recall a generalized 
version of the chain rule for G(u) where G is a Lipchitz function and u a 
Sobolev one. 



(S) 
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3 Mathematical background 



In this section we shall recall some results concerning: the truncature at a 
fixed level and the Stampacchia's estimates. This last technique takes an 
important role in our analysis, moreover we shall need a precise control of the 
estimates. Hence we shall present here the technique with some details and 
developments. As in the rest on the paper we denote by SI C a bounded 
open Lipschitz domain. These properties for f2 are always implicity assumed 
if they are not precised. 



3.1 Truncatures and related properties 

The technique of Stampacchia is based on the use of special test functions 
which are constructed by using some truncatures. We shall recall some basic 
properties of the truncatures used in the paper. An important tool is the 
generalized chain rule (see Theorem [13] in the Annexe). 
Let I > we denote by T\ the truncature function 7} : R — > R defined by 

(I if s > I 

Ti(s) = < s if - I < s < I (22) 
{ -I if s < -I 

Let v G H 1 ^). By applying Theorem[Hwe see that T t (v) G H 1 (O) D L°° (U) . 
Moreover, if we denote by Q v i the set {i£fl s.t. \v (x)| < /} then we have: 



VTj(«) 



Vu in £l Vt i 
elsewhere 



Note that T;(.) truncates both the positive and the negative large values. In 
some cases we need only to truncate in one side. For this reason we introduce 
the semi-truncatures T} )+ and Ti_ defined by: 

T, (t)= l l iH>/ Tr(t)= { ~ l {U< ~ 1 

' \t elsewhere \ t elsewhere 

We then have the decomposition: 7} = o = o 7} +. 
For given V[ G H 1 ^) and s > 0, we shall also consider 

( Vl -s inA+ t 
ip s (v) = vi- T s (yi) = «i + s in (23) 

^ elsewhere 

where we have used the notations: 

A+j := {vi > s}, A~, := {v t < -s} and A s j = A+ t U A~ v (24) 

Let also i/) a ± be the functions defined above (|23|) while replacing T s by T S)+ 
or by T Si _ . It is easy to verify that ip s ,+ (resp. ip s ,+) is in fact the positive 
(resp. the negative) part of ip s . In other words, we have : 

^ St+ (v l ) = ^+(v l ) = (v l -s)l A+ >0, ^ St -(v l ) = ^-(v l ) = (v l + s)l A - <0. 

s,l s,l 

The function tp s has the following properties: 
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Lemma 4 

i) ip s G and = l.A,,,Vwj 

m,) G L°°(dQ) then for all s > \\7v1 \ \ L°°(dn) we have ip s G Hq(Q), 

where 7 : H 1 ^) — > H l / 2 {d^l) denotes the trace function. 
Proof 

Point i) is a direct consequence of Theorem [I3j Property ii) is proved in |14j . 
lemma 3.3 p. 53. □ 



3.2 The Stampacchia estimates 

The Stampacchia estimates is a general method which allows one to obtain 
an L°°-estimate for the weak solution of a large class of elliptic PDEs of the 
second order. The L°°-estimate presented in the original work [26J depends 
on various quantities related to the PDE problem studied, but the exact de- 
pendency is not established. In our analysis we need a precise control of the 
L°°-estimates with respect to some quantities (in particular with respect to 
the diffusion coefficient of the PDEs). Hence in the following we take over 
and detail the technique in order to obtain a more precise L°°-estimates. 

The Stampacchia estimates are established by using the test functions ip s 
(or defined previously, where in this case v (resp. v{) is a weak solution 
of the problem (resp. the sequence of problems) considered. 
For technical reasons we need a classical result concerning some relationship 
between L r functions and linear form on Sobolev spaces: 

Lemma 5 Let 1 < r < 00 and E G L r (£l). Then 

Nr 

E G W-^ip), with (3 = r* = — — -, (25) 
and there exists E G (L^(f2)) 3 such that: 

I Eip= I E-V^ V^P(fi), \\E\\ {L , m3 <\\E\\ Lr{n) . (26) 



Moreover we have: 



N 

r > — (3 > N. (27) 



Proof 

Property (|25[) is easy to prove by using the Sobolev injection Theorem together 



l,p ;f * Np 



N-p - 



with the Holder inequality: 99 — > Eip is a linear form on W ' p if p* : 
^-j- . This last condition holds for p = P : = r ^ N ^_ N ■ Hence P' = pj (p — 1) = 

iVr _ * 
N—r 

We next obtain (f26|) by using a classical result (see [2] Proposition IX. 20). 
Finally, if we assume that r > then /3 = r* > N _^/ 2 = N. □ 

The Stampacchia technique works in two steps: the first one is dependent 
of the problem (or the sequence of problems) studied and the second one is 
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independent of it. Here the purpose is to present the key ingredients of these 
two steps. Because the first one is dependent on the problem studied we 
cannot present it here in its enterity, but we will consider a simple problem 
which contains the main technical points (in fact this introductive presentation 
will be useful to treat a more complicated class of problems in Section 4) . Let 
(vi) C Hq(Q) be a sequence of functions satisfying: 

/ v l Vv l V<p= [ gup, V<p£H%, (28) 
Jn Jq 

where (vi) is a given sequence of strictly positive bounded functions and (gi) C 
L r (Q), with r > y- Let also mi, Mi denote the bounds from above and below 
for i>i, that is: 

< mi < vi < Mi < oo, a.e. in Q. 

-Step 1 

By testing (|28j) with tp s = ifi s (vi), we obtain: 

mi [ | W s | 2 = / g^s by 1 ~ E f EM., (29) 
Jq Jq Jq 

with Ei G lP, \\Ei\\ L fi < \\gi\\^ for some (3 = (3{r) > N. 

Recall also that Supp ip s (vi) C Ai s . Hence by using the Holder inequality we 
obtain: 

f EiV^ s < HVVsM / |^| 2 )^ < ||V^|| L2 ||^|I^(/ i)^ 7 
Jn Ja 1s J A i>s 

< \\Vif> a \\&\\E l \\ LP \A i .\ e sr. 



Consequently (|29j) leads to: 

2 N — 2 

H^llHi (n)2 < C\A hs f, whith C> 0, $ > - = -j^-. (30) 

This is the key estimate needed to pass at the second step which is inde- 
pendent of the problem studied. 

Note that with the particular sequence of problems ()28p chosen here, the 
constants C and are: 



lair $ = t 

2 ' o 

mf p 



C= }MHVL i cJ> = ^_A (31) 



Hence <I> does not depend on I, s. Moreover if we assume that {gi) is uniformly 
bounded in the L r -norm and that (mi) is uniformly bounded from above by 
a strictly positive constant, then C is also independent on /, s. 
This is an important point because as we will see it hereafter, an estimate 
(I30D with C and <I> independent on I, s leads to a uniform L°° bound for (vi). 
-Step 2 

Assume that we have obtained ()30p . We can obtain an L°°-estimate for vi as 
follows. 

Let 2* = -J^2 De the Sobolev exponent associated to 2 in dimension N . By 
using the Poincare-Sobolev inequality we have: 

(/ \H r ) 2/T < Ci\m 2 H i in) i Ci = CidQI) (32) 



n 



Let now t > s. It is clear that A\j C A^ s and consequently 

(/ \H 2 7 /r >(f \H 2 7 /r >(f \t-sff' T >\t-s\ 2 \A ut f T 

■JA,s JA Kt JA lit 

(33) 

We set 

X(s):=\\ s \, Vs>0. (34) 

For fixed Z, x is a decreasing function, and from the estimates (}3"0)) - (f3"3"j) . we 
obtain 

x(t)<C 2 \ X (s)\ 2 **/ 2 (t-s)- 2 * Vt> S >0, C 2 = (C 1 Cf/ 2 . 

Recall that we have assumed in d3D} that $ > Jjr. Hence 2*<I>/2 > 1 and by 
using Lemma 4.1 in |26| we obtain: 

x (d)=0, d = 2^2C 2 1/2 *|0|t"^ < oo. (35) 
This property tells exactly that: 

MU<» ( n)<d, <* = C 3 ($,iV,N)^- (36) 

In particular d does not depend on I if the constants C and <3? appearing 
in (|30p are independent of /. For instance with the particular sequence of 
problems (|28j) the constants C and $ are given by (f3TI) . and if we assume that 
IIs'zIIl 7 ' < K, m| > m > we obtain: 

d = C—, C = C(n,N,r). 
m 



Remark 6 



// you are only interested in obtaining an uniform majoration or mino- 
ration for vi then instead of A30\) it suffices to have 



£l& oW < C+lA+f*, with C± > 0,$± > ^4. 



In fact in this case we consider x±( s ) := I^^J instead of f^p . This 
function is decreasing and we obtain again t35\). But now this property 
tells: 

±v t (x) < d ± a.e. inn, d ± = C ± (^> ± ,N, \0,\)VC±. 

Let again (vi) be a sequence of functions satisfying A28\) and assume that 
mi > m > 0. Then we have: 

gi < hi and H := sup H/i^lli/ < oo =^ vi < d\ a.e. in O, (38) 
i 

gi > ki and K := sup ||&;||£ r < oo vi > —d 2 a.e. in fi, (39) 
l 

with di,d 2 > 0, dt = C-, d 2 = C-, C = C($, N, |0|). 

mm 
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The proof of 138\) and \39}) are obtained by taking over the first step of 
the technique of Stampacchia, but by using ijjf- as test function instead 
ofip s . 

In fact, the function ipf is positive. Hence instead of h29\) we have: 



m / |V^T = / 9i^J < / W = / E lt+ Vift, 
Jn Jn Jn Jn 

with Ei + G L@, \\Ei^ + \\ L p < ||/i;||l<' < H. Consequently, in this case 
we obtain for the function tpf , with C + = and the uniform 
majoration \38\) follows. 

The relation \39}) can be proven by using as test function in 128\) . 
In fact we remark now that ipj is negative. Hence instead of \29\) we 
obtain: 

m [ IVCI 2 = / 9li>7 < I W = / Ei,-V4>7> 
Jn Jn Jn Jn 

with Ei G It, \\Ei < ||A;j||i,r < K. Consequently we now obtain 
JffTy for the function ip~ , with C~ = This implies: —vi < d2 and 
consequently vi > —di- 



4 Approximate sequence and estimates 

Let g G ^^(dtt), we denote by Kg its harmonic lifting, that is: 

Kg G H l (Sl),Kg = g on dU and ARg = on fi. 

We define the functions 6q and cpo by the formula: 

6 := Ka 99 := Kb. (40) 

Hence, by the using the maximum principle (see [2] p. 189 and [11]) together 
with the condition (11511 we obtain: 



< 6 < O < \\a\\ L oo {dn) , 5 < Lpo < ||6||i,oo( a n)- ( 4i ) 
Let now n > 0, (6 n , cp n ) be given and 

cjp\.) :=C e (.,e n (.),v n (.)), C$\.) :=C v (.,e n (.),M-))- 

In order to construct an approximate solution (0 n +i, <p n +i) f° r problem (Q), 
we introduce the following system: 



(Qr 



pu ■ V6 n+1 - div ((y + ^r)V0„+i) = 9^ in 0, 
pu • V^ n+ i - div {{y + g— f— - )Vip n +i) = 9™ in Q, 



G n +\ = a and (p n +\ = b on <9f2, 
where we used the notations: 

gf ] := C 5 - C 3 F9 2 n+1 + C A 6 n+1 D, g^ := -^(Cg^i + C 6 9 n+1 F - C 7 D). 
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For n G N we denote by (H n ) the following condition: 
(Hn) 

Let (/3 max be a fixed real number such that: 



, <^n > 0. 



V?max > ||&IU°°(an)- ( 42 ) 

We shall also consider the condition: (K n ) := (H n ) + (</>„ < <^ m ax)- 

Note that (|4f p shows that the condition (K n ) is statisfied for n = 0. We 
will prove in the sequel that under condition (K n ) we can obtain a weak 
solution (fln+i, tpn+i) f° r problem (Q n ), with moreover (6 n+ i,ip n+ i) satisfying 
the condition (K n+ \). This last property ensures the right definition of an 
approximate sequence. More precisely, we have: 

Proposition 7 Let n G N be given and assume that (Aq) is satisfied. Let also 
(On, Pn) be given and satisfying condition (K n ). There exists r > depending 
only on DATA such that if \\D + \\i J r < r then problem (Q n ) admits at least 
one weak solution (9 n+ i,(p n+ i) C x L°°($l)) 2 . 

Moreover (6 n +i,<p n +i) satifies condition (K n+ \) and the estimates: 

< m i n < 6 n+ i < 9 ma x, (43) 
< (Pmin < fn+l < fmax, (44) 

where ip ma x was fixed in (JW an d fmini &min, G max are positive numbers de- 
pending on DATA, but not on n. 

Remark 8 Proposition is the key result that will be used later on to prove 
Theorem QJ whereas for Theorem [H we shall establish and use a more simple 
version of this proposition (see Subsection \5.2\) . 

In order to prove the proposition we establish intermediate results. 



4.1 Auxiliary results 

Let n G N and (9 n , tp n ) be given and satisfying (K n ). We want to obtain 
(9 n +i, fn+i) by solving (Q n ) and in order to iterate the algorithm we also 
want that (9 n +i,fn+i) satisfies (K n+ i). 

Remark that the system (Q n ) is composed of two coupled scalar elliptic equa- 
tion in divergence form, with a possible singular and degenerate structure. 
Hence the goal of this subsection is to study this last kind of scalar problem. 

In order to do this, we first introduce a weight k : M + — > M + which is 
assumed to be mesurable and satisfying: 

< kq < k < ki a.e. in S7, (45) 

where kq and Ati are two given reals. 

Let also g : £1 x R + — > R be a Caratheodory function (i.e. Vu G M + : 
x — > g(x,u) is measurable, and for a. a. u — > g(x,u) is continuous). 
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Let us consider the scalar problem: 
(S) 



pu-VC-div^+f)VC) =g(x,() inO 
C = Co on dfl, 



where Co G # 1/2 (dO) n L°°(dft), Co > <5 > a.e. in 5Q, is given. We 
always assume that p, u, v, f2 which appear in (S) satisfy their corresponding 
conditions in (Aq). 

Recall that we allow v = in (Aq). Hence problem (S) may degenerate 
(i.e. the viscosity v + ^ may vanish) when C - > oo. Moreover (S) is singular 
(i.e. the viscosity tends to infinity) when C — > 0. 

We want now to find sufficient additional conditions for g that guarantee 
the existence of a bounded positive weak solution for problem (S). Hence we 
shall consider: 

j(x) : = sup \g-(x,u)\e L r (fl), (46) 
ue[o,i] 

g + {x,u) < 7i(x) + 7 2 (x)/t(n), (47) 

where 71,72 G L r (£l) and h : R + —> M + is continuous. In fact, more than 
establishing only the existence of a bounded positive solution for (S), we are 
interested in obtaining some uniform (with respect to K\) bounds from above 
and below and some regularity results. We have: 

Proposition 9 

i) Let k satisfying j^5[ ) and g : £1 x R+ — > R be a Caratheodory function 
satisfying (ffi ), (Jty - There exists a real r > depending on KQ,h such 
that «/||72||L r (f2) ^ T then there exists a weak solution C G H 1 (£l)nL°° (£1) 
for problem (S). Moreover we have: 

-C_ c_ 

^ Cmin ^ C — Cmaxj Cmin — 6 K , Cmax — 6 re ; (48) 



max 



where C depends only on 7, 71, f2, r, iV, Co- particular Cmin o,nd C 
are independent of k\ . 

In addition, the following extended^ maximum principle holds: 

UC[|L~(n) < l|ColU~ ( an)+^^IICIIz°°(n), C = (7(0,iV, r , 7l) . (49) 

22) Assume that in addition u G (L°°(J7)) Ar and Co is Holder continuous. 
Then Q G C 0,a (£l) for some < a < 1. Moreover if <9f2 is of class 
C 2,a g e C ' a (_nxlR+),pu G (C^pj^.KEC 1 ^!!) and Co eC 2 - Q (ffi), 
i/ten C G C 2 ' Q (^) and ^ is a classical solution of (S). 

Before proving Proposition [9] we establish an intermediate result. In a first 
step we consider the change of variable v = InC in (S), and for Z G N we 
introduce a truncated version (Si) of the system obtained: 

, . f puVe Ti ^ -dw((ue Tl ^ + k)Vv) = g(x,e Tl ^) in ft 
{ l) \ u = lnCo on dQ 

We then establish: 



3 when g + = it is a maximum principle 
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Lemma 10 

i) Let k satisfying fl^5| ) and g : £1 x M + — )• M 6e a Caratheodory function 
satisfying f^o] ), f^T) ). T/ien, /or any / G N t/iere eziste a wea/c solution 
v = vi£ H l (n) n L°° /or t/ie problem (Si). 

ii) Let (vi) be the sequence given in i) and I > 1 be a fixed integer. Then 
there exists r\ > such that if the function 72 in |^6| ) satisfies H72IIL'' < 
Ti then we have: 

NU~(o) <c ll7 + 7lllLr , C = C(ZMTA, Co). (50) 

In particular C is independent of n, v and I. 
Proof 



i) By using the divergence formula we obtain, for all w G H l (£l): 

[ puVe Tl( - v) .w = - J Q e T '^ div(puw) + [ e Tl{v) pwu ■ nda 
Jn ' " Jan 



= -fn e Tl(v) Pu ■ Vu; - f Q e Tl ^ div(pu) w. 

=0 

Let vq := \n(lZ£o) and consider the change of variable v := v — vq. Then 
problem (Si) is equivalent to find v G Hq(£1) such that: 

- diva(x,v,Vv) = f(x,v) in £>'(Q), (51) 

where a : x R x -> t w and / : O x M ^ R are defined by: 

a(x,«;,G) = (ve T '( w+Vo) + k(x))G - p(x)u(x)e Ti( ^ o(:c)) , 
/(x,w) = 5 (x,e T ' (u,+t,o(2:)) ). 

We now remark that (|5ip is a quasilinear equation in divergence form. 
Moreover, it is easy to see that / and o~ satisfy the classical growth 
assumptions and a satisfies also the classical coercivity condition. Note 
that: 

(a(x, w, G) - a(x, w, G), G - G) = (ue T ^ w+V0 ^ + k)\G - G| 2 

> k |G-G| 2 . 

Hence a is strictly monotone in the third variable. We then conclude 
(see for instance [1] Theorem 1.5, or [13] Theorem 8.8 page 311) that 
there exists a weak solution v G Hq(Q) for (f5T1) . 

Consequently vi := v + vq is a weak solution for (Si), that is Vu> G Hq(Q): 

[ (ve TlM +k)VviVw+ [ puVe Tl( - Vl) w= [ g(x,e TlM )w. (52) 
Jn Jn Jn 

By applying Theorem 4.2 page 108 in [26] we obtain: vi G -L°°(f2). 
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Let T\ := , > and assume now that ||72||L r < T£- (53) 
rt(e') 

With this additional assumption, we are able to obtain a useful esti- 
mation for \(V\ 1 1 loo .Technically we will detail a method due to Stampac- 
chia (see Subsection 13.21 for the notations and for an introduction of the 
method. Here only the first step of the technique will be developed fur- 
ther). Let s > | In 1 1 Co 1 1 l°° (an) I; we consider the function ip s = vi — T s (vi). 
We have (see Lemma ED ip s £ H^(U) n L°°(Q), and by testing ([52]) with 
ips we obtain: 

f (y^M) +K )|v^| 2 + I pu-Ve TlM ip s = f g(x,e TlM )^ s . (54) 
Jn Jn Jn 

* » ' V V ' 

= :I =:II 

• We will now evaluate the terms I and II. 

The term I is simplified by writing one of its integrand factors, namely 
Ve" r '( s '^ s as a gradient. More precisely we have Ve ri '"''i/) s = V(j, with 
Q E H 1 ^) n L°° (see Lemma HH in the Appendix). Hence by applying 
the divergence formula we see that I vanishes: 

I = I pu- V(i dlv ' = rmula — Q di v p U _|_ / pQ u . n d a = o. (55) 
Jn Jn v — v — ' Jdn 

= S v ' 

=0 by (O 

We next estimate the term II: 

11= [ g(x,e T '^ s < [ g + (x,e T '^)i> s + [ g~ (x,e T '^)^ s . 
Jn Ja+, Ja: 



~-Ih -=Il2 



Remark now that on A s l we have vi < — s < 0, which implies that 
e T d v i) < i. Consequently by using the assumption (|46p we obtain: 



Ih < / Tl>,l A - ■ (56) 
Jn s < 1 

The term II\ is majorated as follows: 



Ih < / (7i + I2h(e l ))^ s < /( 71 + 1)^1 A+ (57) 
by g7J JA+ t by 1531 Jn 



> At this point by using the estimates ([SSD , ([5BD and (|57p together with 
54"|) and (|4"5]) . we obtain: 

/ |VV.| 2 < / (7l A - +(7l + l)V)^ s - (58) 



Note that E G L r (J2) and: ||£|| L r < C , C = \\-f\\ L - + ||7i + l||z,r. 
On the other hand (see Lemma[5|) there exists E G (L^(fi)) 3 satisfying: 
||-E||(L' 9 ) JV ^ ll-^l|L r ; and JqE^P = ^ £ R- ecan a l so that 
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we have assumed in ([8]) that r > — which implies (3 > N. By again 



using the Holder inequality we obtain for (p £ Hq- 



P-2 

EVif < \\E\\( L gjN\\ip\\ H i\Swpp ip\ 2 <3 . 



Consequently (|58j) leads to: 

KolhM^ri < / EV^ S < Col^il^ll^lljyi 

Jn 

Young ineq. Kq C% . A .§=1 

Let $ := > ^r^. We have obtained the estimate: 

ii^ii 2 Hol <c |^r, co=^. 

By now using the Stampacchia estimates (see Subsection 13. 2\\ we obtain 
the existence of a real A independent of I such that \A\i\ = 0. Hence: 
INU~(n) <A, A = C(\\( \\ L ~ (an) ,\n\,N,r)^ 

□ 

Proof of Proposition [9] 

i) Existence and estimates. 

Let (vi)i>i be the sequence given in Lemma [TUl Let also I > 1 be given 
and n := j^rj- We assume that ||72||L r < Tj. It follows from Lemma 
ITUlii) that \\vi\\ L oo < K, where K = K(k , DATA) (K independent of I) 
is the integer defined by K = [A] + 1. 

Let now r := tk and assume that ||72||L r < t. Then we have: ||vk-[|x,oo < 
If. Hence Tk{vk) = u if ■ On the other hand satisfies (Sk), that is: 

puVe VK - div((ue VK + k)Vv k ) = g(x, e VK ) in tt 
vk = In Co on 

Let £ := e° K . We have (see Theorem 1131 in the Appendix) Q £ H 1 ^) n 
L°°(fi) and V£ = QVvk- Consequently ( is a solution of problem (S). 
Moreover ||C|U« (n) < e K = C(k ,DATA). 

On the other hand vk > —K almost everywhere, implies that Q > e - ^- 
a.e. in f2. Hence we obtain ([35]) by setting £ m i n = e~ K and C max = e ■ 
• The estimation (|49p is obtained as follows. 
By using the test function ipf instead of tp s we obtain: 

Cmax<A i: =e G «o . (59) 

This last estimation is only a first step in order to obtain the majoration 
for C maiX announced in (l49l) . 

In fact, let h := v + 5. We have 0<^+-5 2 -=:k <k<^ + 7 s - < 00, 

^ ^min 

and 

pu • VC - div(«;VC) = g~ + g + in 
C = Co on dQ, 
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we can then consider the decomposition C = Ci + C27 where Ci (resp. C2) 
G H x (n) n L°°(n) satisfies the following problem (Si) (resp. (S 2 )) : 



J pu • VCi — div(KVCi) = g in £1 
1 Ci = Co on dQ 

f pu • VC2 - div(KVC2) = g + in 
\ C2 = on dn 

Note that the second member in the PDE in (Si) is negative. Hence, 
by the maximum principle (see [26] p. 80 or [2j p. 191 for a simplified 
situation) we obtain: £1 < ||Co||L°°(9fi) a - e - i n ^- 

By using next the Stampacchia technique (see again Subsection 13.21 
Remark [U]) we major the function £2 as follows: 

,||<7 + ||rr a H 7 lH-L r ||<7 + II rr 

C 2 < C' W " L < e C ^^ W , " L , C = C(ZL4T,4). 

K by (J59J fc o 

This leads to the majoration ([49 p . 

ii) Regularity results. 

• If we assume that u G (L oo (0)) 7V then pu G (L°°(n)) N . Moreover by 
using the estimates (|48p the diffusion coefficient v + ^ is bounded from 
above and below and g(x,() £ L r (12) with r > 3/2 fixed. Hence (see 
Lemma|5|: g(x, Q G with (3 > N. Consequently, by using the 
De Giorgi-Nash Theorem (see pj Th. 8.22) we obtain: C G C°> a {Ti), for 
some a > 0. 

• Assume that in addition we have: 

dQ, Co are of class C 2 ' a , g G C°' a (Ti x M+), 
puG (C^m f, K G C^' a (TT). 

We have proved in the previous point that C G C 0,a (Q). We now iterate 
the Schauder estimates as follows. In a first step we see that g(x, C) and 
v + ^ are in C 0,a (Q), and by applying Theorem 2.7 in [4] we obtain £ G 

C 1,Q (0). Consequently we now obtain (see Appendix B in [8]) g(x, C), v+ 
^ G C 1,a (S7) and by using Theorem 2.8 in [4] we finally obtain Q G 
C 2,a (Q). Hence C is a classical solution of (S). 

□ 



(Si) 
(S 2 ) 



4.2 Proof of Proposition [7] 

Let n G N, and 9 n ,ip n be given. We assume that condition (K n ) is satified. 
Recall that this implies in particular: (p n < <p max , where cp max > ||fr||z,°°(<9n) 
was fixed in (j42j) . Hence, let e := (p max — ||^||L°°(af2) > 0. 
- Step 1: We introduce 

£r(n) 

K (n) . = ^( X; n ) . = Q. _ C3i?u 2 + C4£)u 

<Pn 
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Hence the first subproblem in (Q n ) reads as: 

(Qn-l) 



pu ■ V9 n+1 - div ( (v + ^-)V0 n+1 J = 50(2;, 0n+l) in 



n +i = a on 50 
Moreover, it is easy to verify that: 

0<^< K W<||^|| LOO <oo 

= :ko 

gj(x,u) = C 5 + C A D + u, gg(x,u) = -C 3 Fu 2 + C 4 D~u. 

Note also that Ko is independent of n. Hence we can apply Proposition (9ji) 
(take C = O n+l ,K = K<- n \g = gg,ji = C 5 ,~f 2 = C A D + ,h(t) = i, 7 = C 3 F + 
C4|D~|,Co = a). We obtain the existence of tq > independent of n such 
that if ||-D + ||ir < ro then problem (Q n .l) admits at least one weak solution 
9 n+1 € H 1 ^) n L°°(n). Moreover we have 

< e~ Cipn ^ < 9 n+ i < e c ^" < oo, C = C(DATA). (60) 

- .Step 2: Let now 
/"r(n) 

* (n) := 4 n) (^,^) = = -Vn(C 8 Ci + Ce^n+iF - C 7 D). 



With these notations, the second subproblem in (Qn) reads as: 
(Qn-2) 



pu • V<^„ + i - div Mz/ + ^7)Vv?n+i) = fl$^ in 
(fn+i = b on 90 



We verify that: 



\\r<( n ) II 

0<7 ^<«(«) < l^Jl^l <00 , g+ = <p n C r D+<<p maal C 7 D+, 



'max 
= :kq 



= -MCaOZ+i + Osfln+i-F - C 7 D~), \g-\ < C(DATA, ^ max ). 

Hence we can apply the proposition [SJi) (take now £ = <p n+ i,K = k^ n ',g = 

9(p\ 71 = <AnaxCV.D + , /i = 0, 7 = (5, Co = &)• Then we obtain the existence of 
a weak solution^ 6 ff 1 (0) n L°°(0) for problem (Q n .2). Moreover we 

have 

< e ° < ip n+1 < e° < oo, (61) 



where C,C depend on DATA but not on n. Moreover, by using (|49j) we have: 

lk«+ilU« (n) < ||6|U<- ( en) + C ,, ||^+|Ure c * ta "|^„+i|U-(n), C" = C"(£MTA), 
< l|fc||L^(on) + g1|g 7 ||L^e c ^J |I) + || L H|^ + i||L^(n)- (62) 



4 at this stage there is anymore additional condition needed for ||Z? + ||z,r because h = 
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Assume now that 

ll-t^llz/- < r := min(r , — ). 

'/'max-" 

Then O) leads to 



llVn+lllz^fi) < ll^llL°°(8n) H ll9 3 n+l||L°°(C) ) 

'/'max 

and it follows: 

max • 

V-'max 

- Final Step: if we assume that ||£>+|| L r < r then by using the results es- 
tablished in the previous two steps, we conclude that there exists a solution 
tfn+i) £ {H l n L°°) 2 for problem (Q n ). Moreover this solution satisfies 
(K n+1 ) and ([Ml),® hold. □ 



5 Proofs of the theorems 

We begin by a lemma: 

Lemma 11 Under the assumptions of Proposition^ we can extract a subse- 
quence (still denoted by (9 n ,tp n )) such that 

6 n ^9,tp n ^ip in L°°(0), n ^0,tp n ^ip in H 1 ^), (63) 
' ^-,C { 9 n) -> C e (x,9,<p),CW -»• C v {x,0,<p) ^ L p (0), p < oo (64) 



Proof 

The first properties in (|63p follow directly from Proposition [7J By next using 
n — 0o as test function in (Q n .l) and — ipo as test function in (Q n .2) 
we obtain a uniform bound for (0 n ) and (<p n ) in the iJ 1 -norm. Hence the 
second properties in ([631) follow. Finally property ([64"1) is obtained by using 
the dominated convergence theorem. In fact we have: 



1 1 



, C^ n) -> C e (x, 0, p), -> C^s, 0, a.e in 0, 
^ n ^„ yip 

' ' < t; - , |Cf } | < sup Qfov.u;), |CW|< sup C v (x,v,w), 



Vn^Pn t'min'/'min (u,m)e/C ' (v,w)&K 

where /C = [0, 6» max ] x [0, </? max ]- D 

5.1 Proof of Theorem [I] 

By using (j63[) together with (|64p we obtain: 

pu.V0 n+ i pu.V0 in L 1 (H) 

(u+-^—)V0 n+1 -± (f+/)V0 in q < 2 

0n+l¥>n 0<P 

pu.V<p n +i — pu.Vtp in L 1 (Q) 
(v + tt^ ) Vl Pn+i -± (v + -^)Vip inL*(fi), q<2. 
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Moreover by using ([63]) together with the property 6 n+ \ > 9 m { n > we obtain: 

9 2 n+1 -> 9 2 , 9 n+1 0, 0^ 0" 1 in L p (0), p < oo. 

Hence we can pass to the limit in the approximate problems (Q n ). We obtain 
a weak solution (0, ip) for problem (Q). That is for all ip G V(Q): 

I pu-V9i>+ [ (v + ^)V8-ViP= [ (C 5 -C 3 F8 2 + C 4 8D)iP 
Jn Jn Up Jn 

[ pu ■ Vip V + f(v + ^)Vip • VV = I -piCsV- 1 + C 6 9F - C 7 D)ip 
Jn Jn U P Jn 

= a,{p = b on dO, 
Moreover this solution satisfies: 

9,<p G H 1 ^) n L°°(Q), 9,ip> min(0 min , ^min) > a.e. in Q. 

5.2 Proof of Theorem H 

When N = 2 the function i* 1 has a stronger property of positivity (see lemma 
[12] in the Appendix): 

D 2 

F >^~- ( 65 ) 

We will see that this last property allows one to obtain a weak solution for 
problem (Q) under the assumptions (Aq) and (Ai) but without assuming a 
low compressibility condition of the form (|20p . 

In order to prove this result, we take over the proof of Proposition [7J with 
slight modifications: if {9 n ,Pn) is given and satisfies (H n ) (it is not useful to 
consider (K n ) here) then problem (Q n ) has at least one solution (6 n+ i,p n+ i) 
satisfying in addition (H n+ \) and the estimates (I43p . (l44l) . 

- .Step 1: By using property ([65]) we ma j ore as follows: 

g£(x, u) = C 5 + C A D + u - C 3 Fu 2 < C 5 + C A D + u - ^{D+uf 

< C 5 + D+u(C 4 - -^(D+u) <C 5 + —±. 
^ 1 ' ° 3 

<0 if D+u>?£± 

Hence we have here estimated grj" independtly of the second variable. We 

then apply Proposition [9]i), but now we take 7i = C5 + -^j— instead of C5 
and h(t) = instead of h(t) = t. It follow that there exists (without any 
condition on ||D + ||xr because h = 0) a weak solution 9 n+ \ G i? 1 (0) n L°°(0) 
for problem (Q n .l), with the estimate: 

< e -C||vn|| £ « < Q n+1 < e C\\vn\\ L °° <OQj c= C(DATA). (66) 

- S'tep 2: By taking over the arguments presented in the proof of Proposition [7] 
we see that problem (Q n .2) has at least one positive solution ip n +i G H (O) fl 

i°°(n). 

Hence at this point we have obtained a weak solution (9 n +\, p n +l) f° r (Qn) 
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satisfying in addition {H n+ \). It remains to prove that the estimates (|4"3"]) .([I4" P 
hold. We have made a first step in this direction by proving (|66p . We will 
now prove: 

< e-CK+ilk- < ^ n+1 < \\b\\ Loom , C = C(DATA). (67) 

In fact, by using the additional assumption (Ai) we majore the function 
as follows: 



, x b y E3 Cf. 

g$\x) = -MCsOnli + C 6 6 n+l F - C 7 D) < -^(Cg^ + -f dn+lD ~ ° lD) 

< -^(cs + ^(^ t+ i^) 2 - c 7 (^ +1 D)) = -^V(9 n+1 D), 

Vn+l o Vn+1 

with V(X) := ^X 2 — c-jX + eg. We remark that the discriminant A of is 
negative: A = c| — |c6Cs = —4.864 * 10 -3 < 0. It follows that V is positive 

and consequently gjp is negative. Hence by applying Proposition 0i) with 



now g + = we obtain ([6j 
- Final Step: By using ([66]) together with ([67]) we obtain the estimates (|4"3"j) 
and (|44p . Hence we have recovered the conclusions of Proposition [?J The 
remaindee of the proof for Theorem [3] is exactly the same as for Theorem [TJ 
we can extract a subsequence with the properties ([63j) - ([63|) . These properties 
are sufficient to pass to the limit n — > oo in (Q n ) and Theorem [2] follows. 

5.3 Proof of Theorem [3] 

Let (9, <p) be a weak solution for (Q) in the class S and consider the notations: 

gg(x, u) := C 5 - C 3 Fu 2 + C 4 Du, g v (x, u) := -^fl" 1 + C 6 6F - C 7 D). 

c c 

i) It suffices to remark that the coefficients — & and are bounded from 



above and below, and go, g^ are Caratheodory functions satisfying ([46|) 
and (I47p . Hence we can apply the first point in Proposition [9]-ii) in each 
equation of (Q). We obtain: 6,(p G C°' a (£l), for some a > 0. 

ii) Assume that in addition we have: 

dQ, a, b are of class C 2 ' a , F G C°' a {Tl), 

pu G (C 1 ^))^, C e , C v G C 1 '"^ x (M+) 2 ). 

We remark now that the conditions in the second part of Proposition 
0-ii) are satisfied for each equation of (Q). Hence 0,ip <E C 2,a (£l) and it 
is a classical solution of (Q). 

6 Appendix 

6.1 Derivation of the 9 — if model. 

The model is constructed from the k — e one which takes the form: 

c k 2 k 2 2 

dtk + u • Vfc- —div (0 + p— )Vfc) = c^— F - -/c£> - e, (68) 
p e e 3 

r k 2 2ci f 2 

d t e + u • Ve - —div ((u + p— We) = c x kF -eD - c 2 — , (69) 

p e 2>c u k 
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where D(x,t) := divu(x,i), F(x,t) := ^|Vu+ (Vu) T | 2 - ^D(x,t) 2 > (see 
the next subsection ) and c u ,c £ ,ci, c 2 are generally taken as positive constants. 
Their usual values are (see [23] p. 122): 

C = 0.09, c £ = 0.07, a = 0.128, c 2 = 1.92. (70) 

We then consider the new variables 

6=*, <p = k a e\ (71) 

with a and (5 to be chosen appriopriately. Let Dt denote the total derivative 
operator. By using ([68]) together with ([69]) we obtain an equation for 9: 

89 lk 
D t 6 = — + u-V9 = -D t k - -zDte = -c 3 9 2 F + c A 9D + c 5 + Diff e , 
at e e l 

2 c i 

c 3 = ci - c v , c 4 = -( 1), c 5 = c 2 - 1, (72) 

where Diffg denotes the collected terms coming from the viscous one in the k 
and £ equations. The equation for ip is obtained in the same way: 

k 2 2 

D t tp = ak a ~ l e p D t k + p^e^Dts = ak a - l e^{c v —F - -kD - e) + Dih% 

s 3 

2n f 2 

+ pkfeP-^cikF -eD - c 2 — ) + Diff2^ = Fk a+1 e p ~ l {ac v + fa) 

- k a e?D-(a + P-) - k a - l e^ l (a + f3c 2 ) + Diff,,, 
3 c,. 



= -v(c 6 9F - c 7 D + C$9- 1 ) + Diff,,, 

2 ci 

C6 = -ac^ - /3ci, C7 = —-(a + /3 — ), cs = a + /3c2, (73) 

3 c„ 

where Diff^ = Diffl^ + Diff2 ¥ , is the sum of the terms coming from the viscous 
one in the k and e equations. 

The usual constant values for the parameters 03,04,05 appearing in the 
model are obtained by replacing the values ([70]) in the expressions ([72]) . This 
leads to: 

c 3 = 0.038, c 4 = 0.2815, c 5 = 0.92. 

At this stage it remains to choose appriopriately a and j3 in (|7ip . and to 
model the terms Diff# and Diff^. 

It is schown in [23] , page 67, that a good choice in the incompressible situation 
(i.e. when D = 0) is for instance a = —3, (3 = 2. This leads to the following 
constant values: 

c 6 = 0.014, c 7 = 0.104, c 8 = 0.84, (74) 

and this makes the dynamic stable for the equation in (p. I.e. in the absence 
of the viscous part Diff^, we have: DfCp < 0. 

In the compressible situation the authors suggest in [23], page 125 to consider 
another choice: a = —2— ~ —2.88 and /3 = 2 which makes again the dynamic 
stable. 

Nevertheless in this last situation, the variable <p does not have a clear phys- 
ical meaning (whereas when a = —3, j3 = 2 we have ip = and L := tp^ 1 ! 2 
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represents a length scale of turbulence (see [27])). Moreover, a carefully esti- 
mation schows (see Lemma [12] in the next subsection) that when N = 2 we 
have F > \D 2 . This leads to: 

-<p(c 6 eF - c 7 D + C80- 1 ) < -|(| {ODf - c 7 (6D) + c 8 ) = -^V(0D), 

with V{X) := jX 2 — c 7 X + c$. Hence the choice a = —3, /3 = 2 makes again 
the dynamic stable when N = 2. In fact, in this case cq, c?, cs take the values 
(JHD and the discriminant A of V is A = cf - |c 6 c 8 = -4.864 * 1(T 3 < 0. 
Consequently V(X) > and < in the absence of the viscous terms. 

In consequence we point out that the choice a = — 3, /? = 2 is also in- 
teresting in the compressible situation. We shall make this choice in all the 
situations. Our analysis (see Theorem [1]) shows that this leads to a well posed 
model even when N = 3 under an additional assumption of low compressibil- 
ity of the flow. 

The terms Diffg and Diff^ are modelled (see |16j ) by: 

Diflfe = - div((i/ + cgv t )V9), Diff,, = - div((i/ + c^ t )V6), (75) 
P P 

where vt '■= P c v^- = Pj^ is the turbulent viscosity coming from the equation 
of k, and c#, c v are two new parameters for the model. 

The determination of the parameters eg and can be realized in the same 
way as for the determination of the coefficients arising in the k — e model (see 
|23j). In [22 a constant value for both eg and was numerically tuned from 
a simulation of a Poiseuil flow. However better results are obtained if we allow 
cg,^ to be some positive functions (see [24]). In our analysis we allow the 
coefficients to be of a very general form, in particular cg,^ may depend on 
x, 9 and (p. We only assume that they are Caratheodory functions and that 
they satisfy some positivity and boundedness properties (see (fTT]) - (fT9]) . where 
Cg = peg and C v = pc v ). 



6.2 Positivity of the function F 

In this paragraph we will establish some properties of positivity for the func- 
tion F appearing in the models. 

Let .A/f jv(M) denote the vector space of the N-square matrix with real coeffi- 
cients, equipped with the scalar product: 

A : B = ^ Yl a H h Hi yA = (°y)> B = G M N (R). 



Hence \A\ 2 := \J A : A defines a norm on .Mjv(M). 

For a vector field u : O — > R N we classically defines the gradient Vu and 
the divergence divu (=:D) by: 

dui(x) 



Vu: fl^M N (R), (Vu(x)) tJ - 

. OXi 



25 



Recall that the function F was defined by the formula: 

F(x):=i|Vu+(Vu) T | 2 -^(x) 2 , (76) 

and by an easy calculation we obtain: 

F(x) = (Vu + (Vu)*) : Vu - \d(x) 2 . (77) 

This last expression is sometimes chosen (for instance in |23j) to equivalently 
define F. 

The important fact is that we always have F > but moreover, when N = 2 
the stronger estimate: F > ^D 2 holds. These properties are established in 
the following lemma: 

Lemma 12 The function F satisfies the estimates: 

F = ^ ((dim - d 2 u 2 ) 2 + (d lUl - 8 3 u 3 ) 2 + (d 2 u 2 - d 3 u 3 ) 2 ) 

+ (d 2 m + d\u 2 ) 2 + (d 3 m + d\u 3 ) 2 + (d 3 u 2 + d 2 u 3 ) 2 > 0, when N = 3. 
F = (dim - d 2 u 2 ) 2 + (d 2 m + d lU2 ) 2 + \d 2 > \d 2 , when N = 2. 

Proof 

Let N = 2 or 3, and M := Vu + (Vu)*. Then M;j = djUi + diUj and we 
obtain: 

\M\ 2 = Y, {{IdiUif + J2(d jUi + diUj) 2 ) = 4 ^(c^) 2 + 2 2 + 9i ^') 2 ' 

F = \\M\ 2 - \D 2 = 2 J^u,) 2 - |(£ ftu*) 2 + £ + ftu,) 2 

i i i j>i 



The term ^4 is evaluated separately in the cases N = 2 and N = 3. 

We remark that 2(a 2 + b 2 + c 2 ) - |(a + 6 + c) 2 = |((a-6) 2 + (a-c) 2 + (6-c) 2 ). 

Hence, when iV = 3we have 

4 = | ((5m! - d 2 u 2 ) 2 + (d lUl - d 3 u 3 ) 2 + (0 2 u 2 - d 3 u 3 ) 2 ) , 

and we obtain the expression announced for F. 
In the same way, we remark that 

2(a 2 + b 2 )-\(a + b) 2 = 2 -(a-b) 2 + 2 -(a 2 +b 2 ) 

= 2 -(a-b) 2 + l -(a + b) 2 + ~(a - b) 2 = (a - b) 2 + i(a + 6) 2 . 

Hence, when N = 2 we obtain: 

A = (5itii - d 2 u 2 ) 2 + i + d 2 u 2 ) 2 , 
3* „ ' 

and the expression for F follows. □ 
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6.3 On the low compressibility assumption 

We will show here that without any assumption of low compressibility of the 
form (|20p . problem (Q) may be very hard to analyze when N=3, and singular 
solutions or non existence of weak solution may occur. 

When the dimension equals two we have seen in Theorem [2] that a con- 
dition of low compressibility is not necessary. The reason is related to the 
fact that a stronger property of positivity for F holds in this case, i.e. we 
have: F > When the dimension equals 3 such a property does not hold 
in general. 

In fact, let for instance f2 = B^:i(0, 1) and 

3 

P = ([[{ x i + 4))" 1 , u = (ui,U2,U3) t , Ui=Xi+ 4. 
8=1 

Then a simple calculation gives: D = 3 and F = 0. Hence (Q) reads as: 

' pu • V9 - div (V + i£)V0) = 3C A 9 + C 5 in Q 

< pu-Vip- div ((u + fyVtp) = p(3C 7 - CsO- 1 ) in O 
9 = a,tp = b on d£l 

In this situation the problem becomes hard to analyze. Assume however that 
we have obtained a solution (9, tp) in the class S. Then the equation satisfied 
by 9 is closely related to: 

f - div (77W) = 3C 4 e e + <?5 hi O, 
^ \ = a on dn, 

with bounded from above and below. Hence a contradiction can occur 
because the problem (R) may not have any weak solution (see for instance 
O [19] ) . Note that in the considered example p and u satsify all the conditions 
needed in (Aq), except u-n = on d£l but this is not restrictive for the purpose 
here. In fact we can consider the domain £l\ = B^s (0, 2) which contains O 
and we can extend p, u in f2i in such a way that all the conditions in (Aq) are 
satisfied. Hence we obtain an example within the main situation of the study 
but the evocated problems remain the same. 

6.4 A generalized chain rule 

Let u £ VF 1,p (f2) and G : R — > E be a Lipschitz function. We recall here 
some useful properties of the composed function G(u). In particular we 
schall see that G(u) E W 1,p (£l). Moreover in some situations we also have 
£-G(u) = G'(u)^. 

The main result we have in mind is Theorem [13] which is due to Stam- 
pacchia. In particular we point out that the additional assumption G(0) = 
for the Lipschitz function G is only necessary if f2 is unbounded and p ^ oo, 
or if we want that G(u) has a vanishing trace on dVL when u has it (this last 
situation was in fact the case of interest of Stampacchia). 
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Theorem 13 Let G be a Lipschitz real function, $7 C M. N be a bounded open 
Lipschitz domain, and u G W l,p (VL), with p G [l;oo]. We have: 

i) G{u) G W l *{Q). Moreover if u G W^ P {Q) and G(0) = then G{u) G 

ii ) If G' has a finite number of discontinuity. Then the weak derivatives of 
G(u) are given by the formula: 

-±-G(u) = G'(u)-£- a.e. inn. (78) 

OXi OXi 

Proof 

See the appendix in [26] for the original proof or |12| Theorem 7.8 and [9] 
Theorem 4 for alternative proofs and additional comments. 
We also recall that the formula (178j) may be interpreted in some critical points. 
In fact, let (ti)i=i,.. )Tl denote the points of discontinuity of G' and let E{ := 
{x G ft : u{x) = t{\ be the associated level sets for the function u. Let 
1 < i < n be a fixed integer. If \Ei\ > then the formula (|78p has a priori 
no sense in this last set which is not negligible. Nevertheless it can be shown 
(see [26]) that J^- = on such a set. Hence we interpret the right hand side 



dxi 

of (|78p as zero in the critical set Ei. □ 

We now establish some technical results used in the proof of Lemma [TUJ 
Let vi G i? 1 (f2) be given and consider the function h± G L / 1 oc (]R) defined by: 

h±(y) :=eV(y-T 8 , ± (y)). (79) 

We introduce the functions 

cm) 



9±(t) := / h±(y)dy, teM 
Jo 

Cl,±( x ) := 9±(vi(x)) a.e. x G 0. 



We have: 

Lemma 14 The function Q± has the properties: 

(l,± G H 1 ^) n L°°(n), 
V0,± = VeVMtf. 

Proof 

A simple majoration gives 



|0,±(^)| < / \h±(y)\dy, a.e. in 

JO 



hence Q,± G L°°(n). 

We next remark that g± is a Lipschitz function and its classical derivative is 
given by: 

g' ± {t) = h ± (T l {t))l {mi} Vt^l,-l. 



5 thc derivative G of G takes here the classical sense 
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Hence by using Theorem [13] we obtain 6 H 1 ^) and 

VC/,± = h±(Ti(vi))l { \ Vl \< l} Vvi = h±(vi)l { \ Vl \< l} Vv l 
= e^V^l{M<0^) = Ve T '^'Vf . 

□ 
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